Standing and travelling waves in the shallow-water circular hydraulic jump 
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A wave equation for a time-dependent perturbation about tlie steady shallow-water solution emulates the 
metric an acoustic white hole, even upon the incorporation of nonlinearity in the lowest order. A standing wave 
in the sub-critical region of the flow is stabilised by viscosity, and the resulting time scale for the amplitude 
decay helps in providing a scaling argument for the formation of the hydraulic jump. A standing wave in the 
super-critical region, on the other hand, displays an unstable character, which, although somewhat mitigated 
by viscosity, needs nonlinear effects to be saturated. A travelling wave moving upstream from the sub-critical 
region, destabilises the flow in the vicinity of the jump, for which experimental support has been given. 
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I. INTRODUCTION 

The hydraulic jump is a subject of textbooks in fluid dynamics III S B 0, B • With a jet of water impinging vertically on a 
flat surface, it can be seen that from the point of impingement, the water spreads out radially in a thin film, and beyond a certain 
distance, the height of the fluid layer abruptly increases, exhibiting the circular hydraulic jump |@, 0> [I]- This is a common 
occurrence in a kitchen sink, and somewhat differently this phenomenon is also exhibited as a tidal bore moving upstream in a 
river [;4J, which makes the hydraulic jump an engineering problem of considerable practical interest as well. 

Viscosity of the flowing liquid is the key factor that makes the hydraulic jump possible [13. Hill, and in consequence of this, 
a scaling relation for the hydraulic jump radius has been derived by Bohr et al. |9[[, with strong experimental support coming 
subsequently in favour of this scaling relation through a work by Hansen et al. p2ll . Some later works have also shown that 
surface tension has an involvement with the jump in shallow-layer flows llT3llT3,[l5|] . While dwelling on these matters, it would 
be very pertinent to mention that in a recently published work, Rolley et al. piS] have brought forth an interesting experimental 
study in which the hydraulic jump was seen to form in the flow of liquid helium, below its superfluid transition temperature. The 
radius of the jump and the depth profile of the flow were also not altered appreciably at the superfluid transition, as though the 
fluid was still normal. This feature is probably due to the fact that the flow velocity in this situation is supercritical with respect 
to the critical velocity for superfluidity. 

While it is usual to address all jump-related issues in the stationary regime, in a departure from this practice, explicit time- 
dependence in the shallow-water approach to the circular hydraulic jump has been considered in this paper A time-dependent 
perturbation has been imposed on the steady volumetric flow rate. As an immediate and interesting consequence of this line of 
attack, a linearised equation of motion for the perturbation has been found to have a mathematical form that is very much like 
the metric of an acoustic white hole. This compatibility makes it possible to match the critical properties of the flow with the 
horizon condition of the acoustic white hole, with important consequences leading to the formation of the jump itself. 

As regards the stability of stationary solutions, it has been shown that when the perturbation is treated as a standing wave, 
bounded at two separate spatial points in the sub-critical flow region, viscosity will cause a damping of the amplitude of the 
perturbation. A resultant time scale for viscous damping of the perturbation has been argued to be the same time scale on which 
the flow is dissipatively slowed down by viscosity, and using this time scale, along with the analogy of a white hole, it becomes 
possible to establish a physical explanation for the standard scaling relation of the hydraulic jump radius Ql- A study to this 
effect has been reported before il6[l for the one-dimensional channel flow. 
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In the super-critical region, however, the standing wave does not exhibit similarly stable behaviour. Viscosity goes some 
way in opposing the growth of the perturbation in this region, but this divergent aspect of the perturbation has to be saturated 
ultimately by the conspicuously nonlinear effects in the super-critical region of the flow. This feature is in keeping with what 
has been seen from earlier experiments 1 17]. A very intriguing feature that has emerged on including nonlinearity is that there is 
no alteration of the conditions under which it is possible to construct an acoustic metric for the flow. 

As opposed to standing waves, treating the perturbation as a travelling wave leads to a different conclusion about the stability 
of the flow. In this case it can be seen that the outgoing mode of the travelling wave has a growing amplitude in the super-critical 
region of the flow, but it ultimately decays out in the sub-critical flow region. On the other hand, the incoming mode causes 
great disturbance in the flow as it propagates upstream in the sub-critical region, especially very close to the jump. Experimental 
evidence has been adduced for the latter kind of behaviour 1 18]. 



11. A WAVE EQUATION AND THE METRIC OF AN ACOUSTIC WHITE HOLE 



The flow is described by the Navier-Stokes equation and the equation of continuity, with both of them modified for a shallow 
radial flow confined to a plane ]i9j]. The steady flow variables would be the velocity of the flow, v{r), and the local height of the 
incompressible fluid, h{r), with the former having been obtained in the shallow-water theory by vertically averaging the radial 
component of the velocity, u{r, z), over the height, h{r). The boundary conditions for the averaging are u{r) = at z = 0, and 
du/dz = at z = h{r), under the assumption that while the vertical velocity is much small compared with the radial velocity, 
the vertical variation of the velocity (through the shallow layer of water) is much greater than its radial variation ||^. 

The dynamic generalisation of both v and h, and their corresponding dynamic equations may be written as 



dh Id,,, 

5^ + -^(-M = o 
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dt dr ^ dr /i^ ' 
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respectively, with the usual viscous term in the Navier-Stokes equation having been suitably approximated in Eq. (|2]), following 
the arguments offered by Bohr et al. ]9]. The steady solutions of Eqs. ([T]i and (O above lead to the phenomenon of the hydraulic 
jump. 

To carry out a linearised perturbative analysis, it will be convenient to work with a new variable that is defined as / = rvh. 
The wisdom behind this choice may be discerned from the structure of Eq. ([T]), the steady solution of which can be integrated to 
give rvh = Q/2'k, with the integration constant Q being the steady volumetric flow rate. Therefore, the new variable / can be 
physically identified with the time-dependent volumetric flow rate, and its steady solution will be a constant. 

Solutions of the form ti(r, i) = vo{r) + v'{r,t) and /i(r,t) = ho{r) + h' [r^t) are to be used, in which the subscript indicates 
steady solutions, while the primed quantities are time-dependent perturbations about those steady solutions. Linearising in these 
fluctuating quantities gives the fluctuation of / about its constant steady value, /q — rv^h^ — Q/2'it, as 



/' = r {voh' + hov') 



In terms of /', it is possible to write from Eq. ([T]i, 



which, in conjunction with Eq. (O, gives. 
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A further partial differentiation of Eq. (|5]l with respect to t, followed by extracting the linearised terms in v' and h' from Eq. (|2|i, 
and substitution of these perturbed quantities in terms of /', will deliver 
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It should be instructive to examine Eq. ^ in its inviscid limit, i.e. when v — Q. Some studies lT9[ |20ll2lll have revealed a 
close analogy between the propagation of a wave in a moving fluid and of light in curved space-time. On this particular issue. 
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Schiitzhold and Unruh have shown that gravity waves in a shallow layer of liquid are governed by the same wave equation as for 
a scalar field in curved space-time [19]. For an inviscid, incompressible and irrotational flow, the flow velocity has been defined 
to be the gradient of a scalar potential. Perturbing this potential about its background value, under the restricted condition of the 
background flow height being a constant, leads to an effective metric, in which the velocity of gravity waves replace the speed 
of sound in sonic analogues which closely reflect features seen in general relativistic studies LL9, 2QJ. 

The present analysis is devoted to what is essentially a dissipative system (since it includes viscosity). However, in its inviscid 
limit, this system will deliver the same metric obtained by Schiitzhold and Unruh from their purely inviscid model. It must 
be stressed here that the choice of perturbing the constant background flow rate is paticularly expedient, since conservation of 
matter holds good even in a system that undergoes viscous dissipation. It is to be further noted that the background velocity 
and flow height in this treatment are in general stationary functions of space. Extracting the inviscid terms only from Eq. (|6]l by 
setting = 0, it should be easy to ultimately render these terms into a compact formulation that looks like 1.20.1 

9a {r^dpf) = 0, (7) 

in which the Greek indices run from to 1, with the identification that stands for t and 1 stands for r. On inspecting the terms 
in the left hand side of Eq. (|6]l, the symmetric matrix. 



1 Wo 

vq - gho 



(8) 



can be obtained. It is well known that in terms of the metric gap, the d' Alembertian for a scalar in curved space, is given by |l2i 

Av = iV^g'^^dp^) , (9) 

in which g'^^ is the inverse of the matrix implied by gap- Under the equivalence that f"^ ~ -y— gg"'', and therefore, g = 
det (f , an effective metric can immediately be set down as 



1 Vq 

Vq vl -~ gho 



(10) 



which is entirely identical to the one derived by Schiitzhold and Unruh. By matrix inversion the inverse effective metric, gj^^, 
is easily obtained from Eq. ( fTOl l. For this system, in which the spatial dependence is only on the radial coordinate, the critical 
condition, Vq = gho, is also identified precisely as the horizon condition of either an acoustic black hole or an acoustic white 
hole, depending on the direction of the flow. The critical condition for a flow progressing from the super-critical region to the 
sub-critical region, is analogous to the horizon of an acoustic white hole. 

In presenting the perturbative analysis in this paper, however, viscosity has been included in the governing equations, to finally 
arrive at Eq. The presence of viscosity disrupts the precise symmetry of the inviscid conditions described by Eq. (|7]i. This 
implies that the clear-cut horizon condition, obtained from the inviscid limit, will be affected. However, this effect, for small 
viscosity, as Schiitzhold and Unruh have argued, cannot be too drastic. They treated viscosity as a small adjunct effect on the 
inviscid flow, and concluded that the basic properties of gravity waves will not be affected overmuch. One way or the other, the 
most important feature to emerge from the analogy of a white hole horizon shall remain qualitatively unchanged, namely, that a 
disturbance propagating upstream from the sub-critical flow region (where Vq < gho) cannot penetrate through the horizon into 
the super-critical region of the flow (where Vq > gho)- In keeping with this, Volovik has also pointed out that the jump condition 
can be closely related to the horizon of a white hole (21)], a surface that nothing can penetrate. This property of the flow has a 
crucial bearing on the way the hydraulic jump is formed. 

It has also been intriguing to note that the linearised perturbation dynamics in the shallow-water hydraulic jump is remarkably 
identical to the corresponding dynamics of astrophysical accretion on to an accretor (either a star or any compact object) Il22ll23l 
I24II25I1 . with the only important point of difference being that whereas in the former case a disturbance in the fluid propagates as 
surface gravity waves, in the latter case a disturbance propagates as acoustic waves. In all other important qualitative respects, the 
similarities between what is an incompressible two-dimensional hydraulic jump outflow taking place on the laboratory scale, and 
what is a compressible spherically symmetric accretion inflow happening on the astrophysical scale, are very much unexpected. 
Apropos of this matter, however. Landau and Lifshitz have made a general but very significant observation that "there is a 
remarkable analogy between gas flow and the flow in a gravitational field of an incompressible fluid with a free surface, when 
the depth of the fluid is small..." iQ]. 



III. STANDING WAVES IN THE SUB-CRITICAL REGION AND SCALING ARGUMENTS 



The perturbation is treated as a standing wave, which is constrained to vanish at two separate boundaries. Between these 
two boundaries the flow has to be continuous. Since the jump is a discontinuity in the flow, the boundaries of the standing 
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wave perturbation will have to be chosen on one side of the jump only, although Eq. ^ itself holds true over the entire flow 
space. Therefore, the stability analysis is confined to the sub-critical region of the flow only, where vq < Cg, with Cg being the 
steady value of the speed of surface gravity waves, expressible as Cg = gho. As for the boundaries of the perturbation, one 
of them can be the outer boundary of the steady flow itself, where, by virtue of the boundary condition on the steady flow, the 
perturbation would naturally decay out. The inner boundary on the other hand, may be chosen to be infinitesimally close to the 
jump, which can be conceived of as a standing shock front. Going by the analogy of a white hole, this point will be seen as an 
impenetrable surface by any disturbance in the sub-critical region, and, therefore, the amplitude of the standing wave will decay 
in its neighbourhood. Between these two boundaries (which actually define the entire sub-critical flow region) the flow would 
have entirely lost its laminar character, and it would be most suited to deriving some physical insight about the behaviour of the 
perturbation and the influence of viscosity on it. 

Upon using a solution of the form /'(r, t) = (r) exp(— itjt), a quadratic equation in uj is obtained, and it is given by 
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(11) 



This expression is first multiplied by vqp and, between the two chosen boundaries, an integration by parts is carried out. Imposing 
the requirement that all integrated "surface" terms are to vanish at the two boundaries, will finally give 
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A solution for uj can be found by solving the foregoing quadratic equation. Of immediate interest is the real part of the pertur- 
bation, going as 



"o 



dr 



vopI dr 



(13) 



In the sub-critical region the flow height, ho{r), is very weakly dependent on r, and so di{~iu!) ~ ~v/2h'^. Under inviscid 
conditions, the perturbation will have a constant amplitude in time, but viscosity, with all its dissipative implications, causes the 
perturbation to decay out. This decay of the amplitude of the perturbation will go as cxp(— t^i/2ft,o)i and it will also set up a 
time scale for viscosity to have its dissipative effect, which, to an order-of-magnitude, will be given by tvisc ~ h'o/^- 

The time scale on which viscous drag in the fluid will cause a dissipative slowing down of the flow, will also be given by 
tvisc- The information of an advanced layer of fluid slowing down has to propagate upstream to preserve the smooth continuity 
of the fluid flow. However, this propagation can happen no faster than the speed of surface gravity waves, Cg, and in the region 
where vq > Cg, no information, therefore, can propagate upstream B [Till . So a stream of fluid that has arrived later, moves 
on ahead, unhindered and uninformed, till its speed becomes comparable with the speed of the surface gravity waves, and only 
then does any information about a "barrier" ahead catches up with the fluid. At this stage one may define a dynamic time scale. 



that Vq 



yn 



'•/vq, which is the time scale on which the bulk flow proceeds. Setting 



<dyn with the additional requirement 



T, gives a condition for the "news" of the viscous slowing down finally catching up with the bulk flow itself. The 
continuity equation gives a further constraint that rvoho = Q/Stt, and on using this, a scaling behaviour for the hydrauhc jump 
radius, 



Q 



5/8^-3/8^-1/8^ 



(14) 



is obtained. This is exactly the scaling relation obtained by Bohr et al. a result that was experimentally corroborated later by 
Hansen et al. 1 12]. 

The crucial factor that emerges from the above discussion is that for the formation of the hydraulic jump, the two time scales, 
tvisc and tdyn, would have to match each other closely, when the Froude number, JF, is close to unity. Under these conditions, 
a layer of fluid arriving late is confronted by a barrier formed by a layer of fluid moving ahead with an abrupt slowness. This 
slowly moving layer of fluid flowed past earUer in time, and at far radial distances it has been retarded considerably by viscous 
drag. In this situation there could not be an indefinite accumulation of the fluid, and since continuity of the fluid flow has to be 
preserved, the newly arrived fluid layer slides over the earlier viscosity dragged slowly moving layer of fluid, and what is seen 
is a sudden increase in the height of the fluid layer, a phenomenon that is known as a hydraulic jump. It is conceivable that this 
picture also explains the eddies which form immediately after the jump lfl2ll . 



IV. STANDING WAVES IN THE SUPER-CRITICAL REGION AND THE ROLE OF NONLINEARITY 



So far the mathematical treatment was concerned primarily with developing and explaining a linearised perturbative effect 
on a steady background solution of the flow. However, what is actually being studied is a fluid dynamical system after all, and 
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to derive a complete understanding of this flow system, at one stage or the other, it should be necessary to address the issue of 
nonlinearity. It has been established through some earlier experimental and theoretical works iflTi l26i] that the stability of the 
flow in the super-critical region comes under threat, and in this connection it is generally recognised that in one way or the other 
nonlinearity has an involvement. 

In this context the question that has to be first posed is about the possibility of carrying out a self-consistent perturbative 
analysis that encompasses nonlinearity. Going back to the earlier prescription of v{r,t) — VQ{r) + v'{r,t) and h{r,t) — 
ho{r) + h'{r, t), but retaining the nonlinear term in the definition of /', it is easy to show that /' = r {v^h' + h^v' + v'h'), while 
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The two foregoing results are mathematically exact, and on comparing them with Eqs. (|4|i and (|5]l, they are seen to be rather 
intriguing too. They show that regardless of the nature of the background flow — stationary or otherwise — the mathematical 
interdependence among v', p' and /' retain an invariant form. This leads, as it shall be presently demonstrated, to some interest- 
ing consequences so far as nonlinearity is concerned. Following the principle behind the derivation of Eq. (|6]l, it can be shown 
that on retaining the lowest order of nonlinearity it is possible to arrive at a similar relation that is of the form 
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Vq ho 
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(17) 



in which £o is the operator implied by Eqs. (|7]i and (O, and A is a "switch" parameter which can assume values of and 1 only, 
depending, respectively, on whether one does not or does want to retain nonlinearity in the lowest order The most interesting 
aspect of Eq. ( [TtI i. however, is to be seen in £i, which is an operator that derives solely from accounting for lowest-order 
nonlinearity (i.e. with fluctuations upto the second order only). The structure of £i is exactly similar to the form of C^, as 
Eq. Q gives it, with the only difference being that the corresponding metric elements, referred to in this case as f will involve 
time-dependent first-order fluctuations, unlike the case of the stationary background state, given by Eq. ([8|l. Expressed fully 
these metric elements will read as 



fQ/3 



-Wo 



h voh — v' 
voh — v' Vgh — 2vov' 



(18) 



in which h = h' /ho, and the Greek indices, a and /?, both run from to 1, as before. The fundamentally symmetric form of 
can only mean that even with nonlinearity incorporated in the treatment, in the inviscid limit (i.e. for v — 0), it should be 
very much possible to construct an analogue acoustic black hole metric (or white hole metric, depending on the direction of the 
flow), following the mathematical procedure adopted earlier in Section |II] Going by the form of Eqs. ( fTSl l and ( fT&t , it is also 
worth surmising at this point that this invariance will not be disturbed even when nonlinearity upto any higher order might be 
accounted for 

Another aspect of Eq. (fTTT i is that apart from nonlinearity it also incorporates viscosity. The question is to what extent viscosity 
is important in determining the fate of a standing wave perturbation in the super-critical region. Going back to Eq. ( fT2] i. but in 
this instance imposing the requirement that the standing wave will be bounded within the super-critical region only, close to the 
jump itself and at a point very close to the origin of the radial flow, it will be easy to see that for vq > Cg, there will be a growing 
mode for the perturbation. A condition for this instability can be set down as 



vq (wg - cl) 
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(19) 



whose form implies that if viscosity, quantified by the value of ly, is large enough, then the instability could be overcome. To 
this extent viscosity might be seen to oppose instability in the super-critical region. However, this conclusion is derived solely 
within the linearised framework, and what is more, in the super-critical region the flow is largely laminar till the jump position is 
reached. Therefore, the true physical role of viscosity may not be of much significance in this region, which is dominated much 
more by nonlinearity. In that event, so far as the standing wave in the super-critical region is concerned, it would be reasonable to 
drop viscosity (by setting v = 0)m Eq. (fTTI l. and retain only the time-dependent nonlinear effects on the background linearised 
stationary state. 
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Having done so, employing the same mathematical methods and physical arguments that led to Eq. ( fTZt . and making use of 
the fact that f ''^ = f a quadratic dispersion relation in u can be derived from Eqs. ( fTTj i and ( fTSl l. which will read as 



dr = 0. 



(20) 

It will be easy to conclude from this expression that the time-dependent part of the perturbation, given by exp(— iwt), will have 
no complex part and, therefore, there will be a divergent mode of the perturbation arising solely due to nonlinearity. 

While all this can be deduced about the amplitude-dependent phase (a regular feature where nonlinearity is concerned), it 
does not indicate anything precise about the behaviour of the amplitude of the perturbation itself. To do so it shall be necessary 
to adopt an approach that is quite common in studying amplitude equations IuTI l28ll . First the expression for /'(r, t) is to be 
modified slightly and written as /'(r, t) = (,{t)puj{r) cxp(— icjt), with (_{t) being a slowly-varying function of time, containing 
information on the nonlinear time-dependent growth of the amplitude. The behaviour of ^(i) is to be studied at a fixed value of 
r very close to the jump itself, where, as experimental evidence shows, the amplitude is saturated by nonlinear effects iITtIi . The 
justification for imposing the requirement of ^(t) being a slowly-varying function of time lies in the fact that the time scale for 
the amplitude amplification of the perturbation in the vicinity of the jump is much greater than the dynamic time scale for the 
formation of the jump, rj/cg. 

Now Eqs. (flSl l and ( fTSl l connect h' and v', respectively, to /'. Considering the modified expression of /', it should therefore 
be easy to appreciate that both h' ^ (^{t) and v' ~ £,{t). By this argument, since uj also contains h' and v', it should be noted that 
uj ^ ^{t). Following this understanding, using the modified form of f'(r, t) in Eq. (fTTj i, and neglecting all second derivatives 
of ^(t), it can now be argued that the general form of the variation of ^{t) with respect to t will be given by ^ = + b^^, 
with any even power of ^{t) being left out to preserve time-reversal symmetry, when ^ — > The coefficients a and b will 
determine the growth of ^, and with a > there will be an early growth mode, but this growth will be saturated if < 0. That 
this is precisely what happens inside the super-critical region close to the jump, is very much in evidence through experimental 
studies iUtIi . With regard to this it would also be quite worthwhile to draw attention to the fact that for a shallow-layer flow of 
liquid helium, below its superfluid transition temperature, a very striking feature is the occurrence of a standing capillary wave 
(stationary ripples) between the impact point of the jet and the jump radius 1,15,1 . 



V. TRAVELLING WAVES 



In treating the perturbation as radially propagating waves of very high frequency, a condition that is needed to be satisfied is 
that the frequency, uj, of the travelling waves should be much greater than any characteristic frequency in the system. Since ivisc 
appears as a natural time scale in the system, it would be logical to use this time scale and impose the condition that cj ^ ^7isc- 
The spatial part of the perturbation is prescribed as Puj{r) ~ e/, with s itself being given by the power series. 



s{r) 



(21) 



The complete perturbation will be given as /'(r, t) = exp(s — \u)t). Following some algebra, Eq. ( fTTT l can be written in a slightly 
altered form as 
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(22) 



in which Puj{r) is substituted by making use of the power series given by Eq. (l2Tl i. After this the three successive highest- 
order terms in w will be obtained as uj'^, uj and w", all of whose coefficients are collected and summed up separately, and each 
individual sum of coefficients is then set to zero. Solutions resulting from expressions involving uj"^ and uj, are given as 



and 



ko 



k-i = 



In (vQCg) ± ly 



Vo ± Cg 



It 



2(J^± 1) 



dr, 



(23) 



(24) 



respectively, with being the Froude number, — uo/cg. 

It will be necessary to show that all successive terms of s(r) will self-consistently follow the condition (jj~"|fc„(r)| ^ 
uj^('^+^)\kn^i(r)\, i.e. the power series given by Pi^{r) will converge very rapidly with increasing n, so that the series can be 
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truncated after the first few terms. This whole approach of studying high-frequency travelling waves, has also been successfully 
applied in astrophysical flows ll22ll25ll . The convergence criterion can be shown to be very much true for the radial shallow- 
water flow here, considering the behaviour of the flow in both the super-critical (where vq ^ constant and ho ~ r~^) and the 
sub-critical (where vq ~ r^^ and ho ~ constant) regions. In the inviscid limit the first three terms in fc„(r) will behave as 
k-i ^ r, ko ^ Inr and ki ~ r^^. With the inclusion of viscosity as a physical effect, it can be seen from Eqs. ( |23] | and ( l24b . 
respectively, that while k-i remains unaffected, fco acquires a ^/-dependent term that goes asymptotically as r. Similarly, two 
such viscosity-dependent terms, going as r and Inr, respectively, are added to fci. However, since ly has a relatively small 
value for water, and since the frequency of the travelling waves is at the same time very great (implied by w ^ ^visc)' 
self-consistency requirement still holds. Therefore, it should be quite sufficient to truncate the power series expansion of 5(7-) 
after considering the two leading terms only, and leave out fci, which will anyway contibute only to the phase of the perturbation 
and not its amplitude. An expression for the perturbation will, therefore, be obtained as 



which should be seen as a linear superposition of two waves with arbitrary constants yl+ and A_ . Both these two waves move 
with a velocity, Cg, relative to the fluid, one against the bulk flow and the other along with it, while the bulk flow itself has 
a velocity, vq. All questions pertaining to the growth or decay in the amplitude of the perturbation will be crucially decided 
by the real terms derived from kg. For the choice of the upper sign, i.e. for outgoing waves, the amplitude will be amplified 
as the wave propagates for a certain distance outwards, and then it will start decaying. A quantitative estimate can be made 
of the asymptotic behaviour of the perturbation. For the highly super-critical region of the flow, much before the jump, vq is 
approximately constant, while Hq ^ r^^. The inviscid term in Eq. ( l24l l will thus contribute a growth behaviour going as r^/^. 
The growth of the viscous term on the other hand will vary as exp [(r/rsi)'^'^], with 7-si being a suitable scale factor At the 
opposite end, when the flow is in the highly sub-critical region, Hq will be nearly constant, while vq ^ r^^. In this region 
the contribution of the inviscid term to the perturbation will be in the form of a growth going as r^/^, which, however, will 
be suppressed by an exponential decay contributed by the viscous term, going approximately as exp(— r/rs2), with being 
another scale factor Interestingly enough, in the vicinity of the jump, the scale factor, rs2, reproduces the exact scaling form 
of the jump radius derived by Bohr et al. |9]. Quite evidently, the perturbation starts decaying after crossing the region where 
JF = 0.5. However, even when the outgoing wave grows in amplitude, the growth is not unrestrained anywhere in the region of 
physical interest, and although it might disturb the jump itself somewhat while passing through it, the stability of the flow as well 
as the jump will not be too much adversely affected overall. So all that happens in the case of the outgoing wave is that starting 
from the super-critical region, it is completely transmitted through the jump discontinuity, continuing to grow in amplitude for 
some distance even beyond the jump, but decaying out eventually at very large distances. 

This relatively placid state of affairs receives a severe jolt with the choice of the lower sign in Eq. ( l24l i. i.e. for incoming 
waves. It is evident from Eq. ( |24] | that for an incoming wave progressing upstream from the highly sub-critical flow region to 
the vicinity of the jump (where the Froude number, JF, will be close to unity), the viscous term will cause a large exponential 
divergence in the amplitude of the perturbation, and the flow will therefore suffer noticeably from large fluctuations. On the other 
hand, in the super-critical region, it is because of the viscous term once again that the flow will be pronouncedly stable under 
any perturbation propagating upstream. Physically speaking this is what it should be anyway, because in the super-critical region 
the bulk flow proceeds much faster against the propagation of the perturbation, and since the super-critical region itself has only 
a finite spatial extent, the perturbation is not allowed to remain here for a span of time long enough to destabihse the flow in 
any way. Indeed, in the super-critical region just before the jump, the amplitude of the perturbation decays to vanishingly small 
values. This can only imply that any incoming perturbation propagating upstream from the sub-critical region is not allowed to 
be transmitted through the jump discontinuity into the super-critical region, and this is entirely in agreement with the white hole 
analogy presented at the end of SectionHJ 

This theoretical contention has received much support from an experimental work carried out by Kate et al. ifTsIl to study 
hydraulic jumps formed due to adjacent normal impinging liquid jets. In a series of clearly illustrative photographs, these 
authors have recorded that when one of the two adjacent liquid jets is gradually moved along a straight line towards the other 
one, there will come a stage when, with the circular jumps formed by the two vertical jets being close enough to each other, 
the liquid trapped between the two jumps will be noticeably pushed up to stand at a much greater height than the rest of the 
flow. In a manner of speaking, a small wall of water — the liquid used in this experimental study — will be formed. It is 
easy to see that the jump formed by the moving jet behaves as a sub-critical disturbance propagating upstream towards the jump 
formed due to the static jet. This disturbance cannot penetrate through the static jump, and so as it approaches the jump region, 
its vertical height increases dramatically, due to accumulation of the liquid. This point has been amply demonstrated by the 
photograph (taken by Kate et al. ifisll ) in Figure [T] However, the disturbance propagating upstream in this experiment, is not 
a purely axisymmetric disturbance, and hence, there is a lack of symmetry in the formation of the wall of water between the 
two jumps. It is conceivable that if the perturbation propagating upstream were constrained to have been perfectly axisymmetric 
in nature, then one would have obtained a manifestly circular wall of water to be caught between the incoming waves and 
the periphery of the circular hydraulic jump. Relaxation for this situation can be achieved through the perturbation expending 




(25) 
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itself in destabilising the jump. That the flow rate (and connected to it by a scaling relation, the jump position itself) is highly 
destabilised by a travelling perturbation incident from the sub-critical side, is a feature that is in quite close correspondence with 
the observation of Hansen et al. [12] that the jump position starts fluctuating beyond a certain critical value of the flow rate. 

The time-averaged energy flux associated with the travelling perturbation can be estimated by first noting that in an annular 
element of fluid, the energy arising due to the perturbation is given by the sum of its kinetic energy. 



?kin = 2^r^ {ho + h') {vo + v'f dr, 



and its internal energy. 



fint = 27rrp 



fho + h'— (hoe) + -/i'^Tj72 (^oe) 
oho 2 oho 



dr, 



(26) 



(27) 



in which p is the constant density of the fluid, and e is the internal energy per unit mass The first-order terms all vanish upon 
time-averaging, leaving an expression for the energy of the perturbation that involves only second-order terms. This reads as 



port 



27rr- 



hov'^ + 2voh'v' + h'^ 



92 



dr. 



(28) 



The height fluctuation, h' , and the velocity fluctuation, v' , could both be expressed in terms of /', by making use of Eqs. (|4|l 
and (O, respectively. Carrying out a time-averaging over the periodic part of the square of /', under the approximation that to a 
leading order s ~ Luk^i, will lead to an expression for the energy flux of the propagating perturbation as 



pert 



^3 



1 92 

l±2T+-—j{hoe] 



exp ±2i/ 



r h"^ 

Cgflo 



It 



2{T±1) 



dr 



(29) 



with a factor of 1/2 deriving from the time-averaging. Once again it can be seen that the incoming wave will cause much 
destabilisation for the energy flux in the region where the value of the Froude number, JF, approaches unity. 



VI. CONCLUDING REMARKS 



Some general comments would well be in order here, considering the broad appeal that the physical problem of the hydraulic 
jump holds across various apparently unrelated disciplines. In this paper a physical explanation for the hydraulic jump has been 
furnished with the help of the analogy of an acoustic white hole, and by applying standard mathematical tools from stability 
studies in astrophysical accretion. In doing so the most salient fact that has been unearthed is that the invariance required 
to invoke the analogy of acoustic geometry holds good even when nonlinearity is taken into account. So this whole study 
has amounted to a very striking convergence of fluid dynamics (both laboratory and astrophysical), analogue gravity and the 
mathematics of nonlinear phenomena. 

The reverse aspect of this condition is also true. Recent literature has found mention of astrophysical problems being addressed 
in terms of the analogy of a hydraulic jump. The radial profiles of large gaseous baryonic structures, such as galaxies and galaxy 
clusters, follow one power law in the inner part, and a different one in the outer part, beyond a certain characteristic radius. 
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This transition has been seen to be sharp. It does not necessarily imply that the two distinctly different profiles are premised 
on different kinds of physical properties. Rather, very much like the hydraulic jump, the governing equations will have two 
solutions, and both will be availed of simultaneously [29]. Another related astrophysical connection pertains to spiral shocks in 
protoplanetary discs. Spiral shocks, for most protoplanetary discs, create a loss of vertical equilibrium in the post-shock region, 
and this results in a rapid expansion of the gas in a direction that is perpendicular to the plane of the disc. This expansion has 
characteristics that are similar to hydraulic jumps ifs^ . 

While these facts are by now well established, the present work may also broach a related issue. In astrophysics, the question of 
merging black holes is presently being much pondered over. Some light could be cast on this question by the kind of experiment 
conducted by Kate et al. 11811, which has provided support for some of the theoretical arguments in this paper. The collision 
of two adjacent circular hydraulic jumps, both modelled as acoustic white holes, could provide a laboratory test case for the 
cognate situation of merging astrophysical black holes. 
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